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Abstract. We study the structural properties and the collective microscopic dynamics of 
, atoms in the amorphous metallic alloy Ni^^Zr^r at the temperature T — 300 A" by molecular 

dynamics simulations. The calculated equilibrium structural and dynamical characteristics are 
compared with the experimental data on neutron diffraction and on inelastic X-ray scattering. 
. We present the interpretation of observed structural relaxation of the microscopic density 

t-H ■ fluctuations of particles for amorphous metallic alloy in the framework of the recurrent relation 

Q , approach. The results of theoretical calculations of the intensity of scattering I{k, cj) for 

O ■ amorphous Ni^aZrer are in a good agreement with the results of computer simulation as well 

as with the experimental data on inelastic X-ray scattering. 

> 

OO I 1. Introduction 

■ Design and production of new materials are important fields in the continuous development 

I of science and technology. Bulk metallic glasses (BMG) belong to the advanced materials that 

attract the special interest of scientists and technologists because of their unique physicochemical 
and mechanical properties [HO [3]. For example, amorphous alloys containing transition metal 
{Fe, Co, Ni) as a main component have high tensile strength, which exceeds by more than twice 
of their crystalline counterparts. Some of these amorphous alloys also have a high corrosion 
resistance, the excellent magnetic and electrical properties. The unique properties of the BMG 



X: 

■ are related with the local structural disorder in the atomic components of the mixture [4j . 

Amorphous metal alloys Ni^ZrioQ-^ are characterized by the following interesting features [5j: 
(i) the structure of the alloys has strong short-range topological order, which significantly affects 
on the dispersion curve in the amorphous state; (ii) the alloys have a large coherent neutron 
scattering cross section, which is an important feature at the comparison of the simulation results 
and the theoretical calculations with the experimental data on neutron diffraction. In addition, 
some metal alloys (such as the known, Zr4i,2Tii3,sCui2.5NiiQ,oBe22.5 - Vitreloy 1) form the bulk 
metallic glass phase at the critical cooling rates of ~ 1 K/ sec |6]. These features have stimulated 
many experimental and theoretical studies. So, for example, the influence of the concentration 
dependence of the NixZriQQ^^ alloys (x = 5,10,16.7,33.3,0.5,66.7,83.3,90.0,95.0) on the 
processes of crystallization of the bulk metallic glass was investigated by using the computer 
simulation methods in Ref. The authors demonstrated that the highest degree of 



crystallization of the alloy is related with high nickel concentration in the system (in particular, 
A^i66.7-^''33.3-3'lloy), where crystallization in the alloys is related with increasing of the number 
of quasicrystalline icosahedral clusters. Further, the study of the dynamic properties of the 
Ni^^ZrQj-alloy in amorphous and crystalline phase was performed in Ref. [8j by using the 
inelastic neutron scattering experiment. By comparing the spectra of the dynamic structure 
factor, S{k,uj), in amorphous and crystalline phases, the authors have concluded that the origin 
of the collective excitations in Ni-^^Zr^f metallic glass is closely associated with the optical 
excitations which observed in the crystalline phase. In addition, the three optical modes, which 
exist in the crystalline Ni^^ZrQj-alloy, can also be observed on the short time scales in the 
amorphous phase. However, the recent experimental studies have questioned this conclusion [9]. 

In this paper we investigate the microscopical structural and dynamical properties of 
amorphous Ni^-^Zr^Y metallic alloy at the temperature T = 300K to elucidate the mechanism 
of propagation of collective excitations. 

2. Details of simulation 
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Figure 1. (Color online) Partial radial distribution functions of particles for Ni^^Zr^f metallic 
alloy at the temperatures T from 300K (top curve) to SOOOi^T (bottom curve). 

The molecular dynamics (MD) simulation of the amorphous metallic alloy Ni^j.ZrQi was 
performed at the temperature T = 300i^. The investigated system was consisted of = 
Njsii + Nzr = 10976 particles, where Nj^n = 3535 and Nzr = 7441; the numerical density was 
n = 0.0427A . The interaction of atoms was carried out with the potential of interparticle 
interaction [TD] 
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Here a,l3 £ {Ni, Zr} and ^q,/3, Ca,i3, o,a,i3, ha^p - are the potential parameters, taken from Ref. 
|10] . All simulations are carried out in the isobaric-isothermal ensemble (NpT). To maintain 
the system in the thermal equilibrium, the Berendsen's thermostat and barostat were applied. 
The amorphous alloy at the temperature T = 300K was obtained by a rapid cooling from the 
liquid state at the temperature T = 3000-R' with the rate dT/dt = 10^^ K/sec. The integration 
of the equations of motion of the particles was performed using the velocity Verlet algorithm 
with the timestep 10~^^ sec. 




Figure 2. (Color online) Static structure factor of A^iasZrer-alloy at the temperature T = SOOif. 



3. Simulation results and discussion 

The structural features of metallic alloy were analyzed by means of the partial radial distribution 
function: 
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and the partial static structure factor [11]: 
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Here, k is the wave-number. The total static structure factor was found by means of the spatial 
Fourier transformation of the total radial distribution function [121 : 
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where, c is the concentration of particles, / is the scattering length (in the case of neutron 
scattering). Fig. 1 shows the temperature dependence of the particle distribution functions, 
9a,i3ii^), calculated on the basis of molecular dynamics simulation. The figure demonstrates 
that the Ni^^ZrQj-system is completely transformed into an amorphous state with a cooling: 
the distribution function is characterized a pronounced first peak and a splitting of the second 
one [13]. Fig. 2 shows the partial static structure factor obtained on the basis of computer 
simulation while the comparison of the MD results for S{k) with the experimental data on 
neutron diffraction jl4j is presented in Fig 2(d). It can be seen from the figure that the chosen 
model of interatomic interaction potential [TU] reproduces in detail the structure of amorphous 
alloy. 
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Figure 3. Intensity of scattering for the metallic alloy Ni^^Zr^i at the temperature T = 300K. 
The solid line represents the simulation results, taking into account the detailed balance 
condition (6) and the experimental resolution function (7); (o o o) - the experimental data 
on inelastic X-ray scattering [9]. The dotted line shows the experimental resolution function 
R{k,uj). 

The microscopic dynamics of the amorphous metal alloy was analyzed by means of the 
dynamic structure factor S{k,uj). To compare the simulation results with the experimental 
data, we have calculated the intensity of the inelastic X-ray scattering, I{k,uj), which is related 
with the dynamic structure factor by the following relations [15\ : 

I{k,uj) = E{k) j R{k,uj-uj')Sq{k,u;')du;', (5) 

Sq{k,uj) = ^ ^^'^^^^^ S{k,uj). (6) 

Here Sq{k,uj) is the quantum dynamic structure factor, /3 = l/^ksT) is the inverse temperature, 
E{k) is a normalization factor, and R{k,uj) is the experimental resolution function 




Figure 4. Intensity of scattering for metallic alloy Ni^^Zr^-j at the temperature T = 300/'C. 
The solid line represents the results of theoretical calculations, taking into account the detailed 
balance condition (6) and the experimental resolution function (7); (o o o) - the experimental 
data on inelastic X-ray scattering [9]. The dotted line shows the experimental resolution function 
R{k,uj). 



which satisfies the normalization condition 



R{k,uj)duj = 1. 



(8) 



Comparison of the simulation results for the scattering intensity I{k, u) with the experimental 
data on inelastic X-ray scattering is given in Fig. 3. As can be seen from the figure, the 
simulation results reproduce accurately the high frequency features of the X-ray scattering 
intensity in a wide range of wave numbers: they predict the positions, heights, and decrease 
of the side peaks. However, the simulation results do not reveal the experimentally-observed 
low-frequency excitations in I{k,io) at the values of fc = 9 H- 11.4 nm~^. 

Accordingly to the theoretical approach developed [161 IIZ] within Lee's recurrent relation 
formalism |18[ [TU] the dynamic structure factor can be taken in the form 
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Here, S{k) is the static structure factor; Ai{k), A2{k), A3(/c) and A4{k) are the frequency 
relaxation parameters: 
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First two parameters were calculated from the molecular dynamics simulation data and the 
parameters A3(A:) and A4(/c) were taken as fitting values. Fig. 4 shows the results of theoretical 
calculations for the scattering intensity I{k,uj) [Eqs. (5) and (9)] with the detailed balance 
condition (6) and the experimental resolution function (7) in comparison with the experimental 
data on inelastic X-ray scattering [9]. It is seen that Eq. (9) yields a good agreement with the 
experimental data. The results of theoretical calculations reproduce the high-frequency features 
of the experimental spectra on inelastic X-ray scattering for all values of wave numbers: the 
correct predictions for the position as well as for the height and damping of the lateral peaks. 

4. Conclusion 

This paper presents the results of the atomic/molecular dynamics simulations of the Ni^^Zr^j 
bulk metallic glass at the temperature T = 300K. The results indicate on the possibility to 
perform the correct description of the structural and dynamic properties of the amorphous 
metallic alloys containing the polyvalent (transition) metals in the framework of pair-wise 
effective potentials. The simulation results of structural and dynamic characteristics are in a 
good agreement with the experimental data on neutron diffraction and inelastic X-ray scattering 
data. The theoretical model developed recently to describe the structural relaxation in liquid 
metals reproduces the features of the I{k,uj) spectra for Ni^^Ztq^ metallic glass associated, in 
particular, with the appearance and disappearance of the high-frequency peaks. 
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